GhRCD1 promotes cotton tolerance to cadmium by regulating the GhbHLH12-GhMYB44-GhHMA1 transcriptional cascade.

Heavy metal pollution poses a significant risk to human health and wreaks havoc on agricultural productivity. Phytoremediation, a plant-based, environmentally benign, and cost-effective method, is employed to remove heavy metals from contaminated soil, particularly in agricultural or heavy metal-sensitive lands. However, the phytoremediation capacity of various plant species and germplasm resources display significant genetic diversity, and the mechanisms underlying these differences remain hitherto obscure. Given its potential benefits, genetic improvement of plants is essential for enhancing their uptake of heavy metals, tolerance to harmful levels, as well as overall growth and development in contaminated soil. In this study, we uncover a molecular cascade that regulates cadmium (Cd2+ ) tolerance in cotton, involving GhRCD1, GhbHLH12, GhMYB44, and GhHMA1. We identified a Cd2+ -sensitive cotton T-DNA insertion mutant with disrupted GhRCD1 expression. Genetic knockout of GhRCD1 by CRISPR/Cas9 technology resulted in reduced Cd2+ tolerance in cotton seedlings, while GhRCD1 overexpression enhanced Cd2+ tolerance. Through molecular interaction studies, we demonstrated that, in response to Cd2+ presence, GhRCD1 directly interacts with GhbHLH12. This interaction activates GhMYB44, which subsequently activates a heavy metal transporter, GhHMA1, by directly binding to a G-box cis-element in its promoter. These findings provide critical insights into a novel GhRCD1-GhbHLH12-GhMYB44-GhHMA1 regulatory module responsible for Cd2+ tolerance in cotton. Furthermore, our study paves the way for the development of elite Cd2+ -tolerant cultivars by elucidating the molecular mechanisms governing the genetic control of Cd2+ tolerance in cotton.

A novel single nucleotide mutation of TFL1 alters the plant architecture of Gossypium arboreum through changing the pre-mRNA splicing.

KEY MESSAGE: A single nucleotide mutation from G to A at the 201st position changed the 5' splice site and deleted 31 amino acids in the first exon of GaTFL1. Growth habit is an important agronomic trait that plays a decisive role in the plant architecture and crop yield. Cotton (Gossypium) tends to indeterminate growth, which is unsuitable for the once-over mechanical harvest system. Here, we identified a determinate growth mutant (dt1) in Gossypium arboreum by EMS mutagenesis, in which the main axis was terminated with the shoot apical meristem (SAM) converted into flowers. The map-based cloning of the dt1 locus showed a single nucleotide mutation from G to A at the 201st positions in TERMINAL FLOWER 1 (GaTFL1), which changed the alternative RNA 5' splice site and resulted in 31 amino acids deletion and loss of function of GaTFL1. Comparative transcriptomic RNA-Seq analysis identified many transporters responsible for the phytohormones, auxin, sugar, and flavonoids, which may function downstream of GaTFL1 to involve the plant architecture regulation. These findings indicate a novel alternative splicing mechanism involved in the post-transcriptional modification and TFL1 may function upstream of the auxin and sugar pathways through mediating their transport to determine the SAM fate and coordinate the vegetative and reproductive development from the SAM of the plant, which provides clues for the TFL1 mechanism in plant development regulation and provide research strategies for plant architecture improvement.

Brassinosteroids regulate cotton fiber elongation by modulating very-long-chain fatty acid biosynthesis.

Brassinosteroid (BR), a growth-promoting phytohormone, regulates many plant growth processes including cell development. However, the mechanism by which BR regulates fiber growth is poorly understood. Cotton (Gossypium hirsutum) fibers are an ideal single-cell model in which to study cell elongation due to their length. Here we report that BR controls cotton fiber elongation by modulating very-long-chain fatty acid (VLCFA) biosynthesis. BR deficiency reduces the expression of 3-ketoacyl-CoA synthases (GhKCSs), the rate-limiting enzymes involved in VLCFA biosynthesis, leading to lower saturated VLCFA contents in pagoda1 (pag1) mutant fibers. In vitro ovule culture experiments show that BR acts upstream of VLCFAs. Silencing of BRI1-EMS-SUPPRESOR 1.4 (GhBES1.4), encoding a master transcription factor of the BR signaling pathway, significantly reduces fiber length, whereas GhBES1.4 overexpression produces longer fibers. GhBES1.4 regulates endogenous VLCFA contents and directly binds to BR RESPONSE ELEMENTS (BRREs) in the GhKCS10_At promoter region, which in turn regulates GhKCS10_At expression to increase endogenous VLCFA contents. GhKCS10_At overexpression promotes cotton fiber elongation, whereas GhKCS10_At silencing inhibits cotton fiber growth, supporting a positive regulatory role for GhKCS10_At in fiber elongation. Overall, these results uncover a mechanism of fiber elongation through crosstalk between BR and VLCFAs at the single-cell level.

Identification and Characterization of the ERF Subfamily B3 Group Revealed GhERF13.12 Improves Salt Tolerance in Upland Cotton.

The APETALA2 (AP2)/ethylene response factor plays vital functions in response to environmental stimulus. The ethylene response factor (ERF) subfamily B3 group belongs to the AP2/ERF superfamily and contains a single AP2/ERF domain. Phylogenetic analysis of the ERF subfamily B3 group genes from Arabdiposis thaliana, Gossypium arboreum, Gossypium hirsutum, and Gossypium raimondii made it possible to divide them into three groups and showed that the ERF subfamily B3 group genes are conserved in cotton. Collinearity analysis identified172 orthologous/paralogous gene pairs between G. arboreum and G. hirsutum; 178 between G. hirsutum and G. raimondii; and 1,392 in G. hirsutum. The GhERF subfamily B3 group gene family experienced massive gene family expansion through either segmental or whole genome duplication events, with most genes showing signature compatible with the action of purifying selection during evolution. Most G. hirsutum ERF subfamily B3 group genes are responsive to salt stress. GhERF13.12 transgenic Arabidopsis showed enhanced salt stress tolerance and exhibited regulation of related biochemical parameters and enhanced expression of genes participating in ABA signaling, proline biosynthesis, and ROS scavenging. In addition, the silencing of the GhERF13.12 gene leads to increased sensitivity to salt stress in cotton. These results indicate that the ERF subfamily B3 group had remained conserved during evolution and that GhERF13.12 induces salt stress tolerance in Arabidopsis and cotton.

Gossypium hirsutum Salt Tolerance Is Enhanced by Overexpression of G. arboreum JAZ1.

Gossypium arboreum possesses many favorable traits including robust defense against biotic and abiotic stress although it has been withdrawn from the market because of lower yield and fiber quality compared to G. hirsutum (upland cotton). It is therefore important to explore and utilize the beneficial genes of G. arboretum for G. hirsutum cultivar breeding. Here, the function of G. arboreum JAZ1 in tolerance to salt stress was determined through loss-of-function analysis. GaJAZ1can interact with GaMYC2 to repress expression of downstream genes whose promoters contain a G-box cis element, affecting plant tolerance to salinity stress. The experimental data from NaCl treatments and a 2 year continuous field trial with natural saline-alkaline soil showed that the ectopically overexpressed GaJAZ1 significantly increased salt tolerance in upland cotton compared to the wild type, showing higher growth vigor with taller plants, increased fresh weight, and more bolls, which is due to reprogrammed expression of tolerance-related genes and promotion of root development. High-throughput RNA sequencing of GaJAZ1 transgenic and wild-type plants showed many differentially expressed genes involved in JA signaling and biosynthesis, salt stress-related genes, and hormone-related genes, suggesting that overexpressing GaJAZ1 can reprogram the expression of defense-related genes in G. hirsutum plants to increase tolerance to salt stress. The research provides a foundation to explore and utilize favorable genes from Gossypium species for upland cotton cultivar breeding.

Genetic and evolution analysis of extrafloral nectary in cotton.

Extrafloral nectaries are a defence trait that plays important roles in plant-animal interactions. Gossypium species are characterized by cellular grooves in leaf midribs that secret large amounts of nectar. Here, with a panel of 215 G. arboreum accessions, we compared extrafloral nectaries to nectariless accessions to identify a region of Chr12 that showed strong differentiation and overlapped with signals from GWAS of nectaries. Fine mapping of an F2 population identified GaNEC1, encoding a PB1 domain-containing protein, as a positive regulator of nectary formation. An InDel, encoding a five amino acid deletion, together with a nonsynonymous substitution, was predicted to cause 3D structural changes in GaNEC1 protein that could confer the nectariless phenotype. mRNA-Seq analysis showed that JA-related genes are up-regulated and cell wall-related genes are down-regulated in the nectary. Silencing of GaNEC1 led to a smaller size of foliar nectary phenotype. Metabolomics analysis identified more than 400 metabolites in nectar, including expected saccharides and amino acids. The identification of GaNEC1 helps establish the network regulating nectary formation and nectar secretion, and has implications for understanding the production of secondary metabolites in nectar. Our results will deepen our understanding of plant-mutualism co-evolution and interactions, and will enable utilization of a plant defence trait in cotton breeding efforts.

Extensive intraspecific gene order and gene structural variations in upland cotton cultivars.

Multiple cotton genomes (diploid and tetraploid) have been assembled. However, genomic variations between cultivars of allotetraploid upland cotton (Gossypium hirsutum L.), the most widely planted cotton species in the world, remain unexplored. Here, we use single-molecule long read and Hi-C sequencing technologies to assemble genomes of the two upland cotton cultivars TM-1 and zhongmiansuo24 (ZM24). Comparisons among TM-1 and ZM24 assemblies and the genomes of the diploid ancestors reveal a large amount of genetic variations. Among them, the top three longest structural variations are located on chromosome A08 of the tetraploid upland cotton, which account for ~30% total length of this chromosome. Haplotype analyses of the mapping population derived from these two cultivars and the germplasm panel show suppressed recombination rates in this region. This study provides additional genomic resources for the community, and the identified genetic variations, especially the reduced meiotic recombination on chromosome A08, will help future breeding.

Over-expression of an R2R3 MYB Gene, GhMYB73, increases tolerance to salt stress in transgenic Arabidopsis.

MYB family genes act as important regulators modulating the response to abiotic stress in plants. However, much less is known about MYB proteins in cotton. Here, we found that a cotton MYB gene, GhMYB73, was induced by NaCl and abscisic acid (ABA). Silencing GhMYB73 expression in cotton increased sensitivity to salt stress. The cotyledon greening rate of Arabidopsis thaliana over-expressing GhMYB73 under NaCl or mannitol treatment was significantly enhanced during the seedling germination stage. What's more, several osmotic stress-induced genes, such as AtNHX1, AtSOS3 and AtP5CS1, were more highly induced in the over-expression lines than in wild type under salt treatment, supporting the hypothesis that GhMYB73 contributes to salinity tolerance by improving osmotic stress resistance. Arabidopsis lines over-expressing GhMYB73 had superior germination and cotyledon greening under ABA treatment, and some abiotic stress-induced genes involved in ABA pathways (AtPYL8, AtABF3, AtRD29B and AtABI5), had increased transcription levels under salt-stress conditions in these lines. Furthermore, we found that GhMYB73 physically interacts with GhPYL8 and AtPYL8, suggesting that GhMYB73 regulates ABA signaling during salinity stress response. Taken together, over-expression of GhMYB73 significantly increases tolerance to salt and ABA stress, indicating that it can potentially be used in transgenic technology approaches to improve cotton salt tolerance.

GhWRKY6 Acts as a Negative Regulator in Both Transgenic Arabidopsis and Cotton During Drought and Salt Stress.

Drought and high salinity are key limiting factors for cotton production. Therefore, research is increasingly focused on the underlying stress response mechanisms of cotton. We first identified and cloned a novel gene encoding the 525 amino acids in cotton, namely GhWRKY6. qRT-PCR analysis indicated that GhWRKY6 was induced by NaCl, PEG 6000 and ABA. Analyses of germination rate and root length indicated that overexpression of GhWRKY6 in Arabidopsis resulted in hypersensitivity to ABA, NaCl, and PEG 6000. In contrast, the loss-of-function mutant wrky6 was insensitive and had slightly longer roots than the wild-type did under these treatment conditions. Furthermore, GhWRKY6 overexpression in Arabidopsis modulated salt- and drought-sensitive phenotypes and stomatal aperture by regulating ABA signaling pathways, and reduced plant tolerance to abiotic stress through reactive oxygen species (ROS) enrichment, reduced proline content, and increased electrolytes and malondialdehyde (MDA). The expression levels of a series of ABA-, salt- and drought-related marker genes were altered in overexpression seedlings. Virus-induced gene silencing (VIGS) technology revealed that down-regulation of GhWRKY6 increased salt tolerance in cotton. These results demonstrate that GhWRKY6 is a negative regulator of plant responses to abiotic stress via the ABA signaling pathway.

Genome-Wide Analysis Elucidates the Role of CONSTANS-like Genes in Stress Responses of Cotton.

The CONSTANS (CO)-like gene family has been well studied for its role in the regulation of plant flowering time. However, their role remains poorly understood in cotton. To better understand the possible roles of CO-like in cotton, we performed a comprehensive genome-wide analysis of CO-like genes in cotton. Phylogenetic tree analysis showed that CO-like genes naturally clustered into three groups. Segmental duplication and whole genome duplication (WGD), which occurred before polyploidy, were important contributors to its expansion within the At ("t" indicates tetraploid) and Dt subgenomes, particularly in Group III. Long-terminal repeat retroelements were identified as the main transposable elements accompanying 18 genes. The genotype of GhCOL12_Dt displayed low diversity; it was a candidate involved in domestication. Selection pressure analyses indicated that relaxed purifying selection might have provided the main impetus during the evolution of CO-like genes in upland cotton. In addition, the high expression in the torus and calycle indicated that CO-like genes might affect flowering. The genes from Group II, and those from Group III involved in segmental duplication or WGD, might play important roles in response to drought and salt stress. Overall, this comprehensive genome-wide study of the CO-like gene family would facilitate further detailed studies in cotton.

Genome-Wide Study of YABBY Genes in Upland Cotton and Their Expression Patterns under Different Stresses.

Members of the YABBY gene family, a small plant-specific family of genes, have been proposed to function in specifying abaxial cell fate. Although to date little has been learned about cotton YABBY genes, completion of the cotton genome enables a comprehensive genome-wide analysis of YABBY genes in cotton. Here, a total of 12, 12, and 23 YABBY genes were identified in Gossypium arboreum (2n = 26, A2), G. raimondii (2n = 26, D5), and G. hirsutum (2n = 4x = 52, [AD]t), respectively. Sequence analysis showed that the N-terminal zinc-finger and C-terminal YABBY domains in YABBY proteins are highly conserved among cotton, Arabidopsis, and rice. Eighty-five genes from eight sequenced species naturally clustered into five groups, and the YAB2-like group could be divided into three sub-groups, indicating that YABBYs are highly conserved among the examined species. Orthologs from the At and Dt sub-genomes (where "t" indicates tetraploid) showed good collinearity, indicating that YABBY loci are highly conserved between these two sub-genomes. Whole-genome duplication was the primary cause of upland cotton YABBY gene expansion, segmental duplication played important roles in YABBY gene expansion within the At and Dt sub-genomes, and the YAB5-like group was mainly generated by segmental duplication. The long-terminal repeat retroelements Copia and Gypsy were identified as major transposable elements accompanying the appearance of duplicated YABBY genes, suggesting that transposable element expansion might be involved in gene duplication. Selection pressure analyses using PAML revealed that relaxed purifying selection might be the main impetus during evolution of YABBY genes in the examined species. Furthermore, exon/intron pattern and motif analyses indicated that genes within the same group were significantly conserved between Arabidopsis and cotton. In addition, the expression patterns in different tissues suggest that YABBY proteins may play roles in ovule development because YABBYs are highly expressed in ovules. The expression pattern of YABBY genes showed that approximately half of the YABBYs were down-regulated under different stress treatments. Collectively, our results represent a comprehensive genome-wide study of the YABBY gene family, which should be helpful in further detailed studies on the gene function and evolution of YABBY genes in cotton.

Genome-wide analysis of BES1 genes in Gossypium revealed their evolutionary conserved roles in brassinosteroid signaling.

Brassinosteroids (BRs), which are essential phytohormones for plant growth and development, are important for cotton fiber development. Additionally, BES1 transcription factors are critical for BR signal transduction. However, cotton BES1 family genes have not been comprehensively characterized. In this study, we identified 11 BES1 genes in G. arboreum, 11 in G. raimondii, 16 in G. barbadense, and 22 in G. hirsutum. The BES1 sequences were significantly conserved in the Arabidopsis thaliana, rice, and upland cotton genomes. A total of 94 BES1 genes from 10 different plant species were divided into three clades according to the neighbor-joining and minimum-evolution methods. Moreover, the exon/intron patterns and motif distributions were highly conserved among the A. thaliana and cotton BES1 genes. The collinearity among the orthologs from the At and Dt subgenomes was estimated. Segmental duplications in the At and Dt subgenomes were primarily responsible for the expansion of the cotton BES1 gene family. Of the GhBES1 genes, GhBES1.4_At/Dt exhibited BL-induced expression and was predominantly expressed in fibers. Furthermore, Col-0/mGhBES1.4_At plants produced curled leaves with long and bent petioles. These transgenic plants also exhibited decreased hypocotyl sensitivity to brassinazole and constitutive BR induced/repressed gene expression patterns. The constitutive BR responses of the plants overexpressing mGhBES1.4_At were similar to those of the bes1-D mutant.

Genome-wide analysis of WOX genes in upland cotton and their expression pattern under different stresses.

BACKGROUND: WUSCHEL-related homeobox (WOX) family members play significant roles in plant growth and development, such as in embryo patterning, stem-cell maintenance, and lateral organ formation. The recently published cotton genome sequences allow us to perform comprehensive genome-wide analysis and characterization of WOX genes in cotton. RESULTS: In this study, we identified 21, 20, and 38 WOX genes in Gossypium arboreum (2n = 26, A2), G. raimondii (2n = 26, D5), and G. hirsutum (2n = 4x = 52, (AD)t), respectively. Sequence logos showed that homeobox domains were significantly conserved among the WOX genes in cotton, Arabidopsis, and rice. A total of 168 genes from three typical monocots and six dicots were naturally divided into three clades, which were further classified into nine sub-clades. A good collinearity was observed in the synteny analysis of the orthologs from At and Dt (t represents tetraploid) sub-genomes. Whole genome duplication (WGD) and segmental duplication within At and Dt sub-genomes played significant roles in the expansion of WOX genes, and segmental duplication mainly generated the WUS clade. Copia and Gypsy were the two major types of transposable elements distributed upstream or downstream of WOX genes. Furthermore, through comparison, we found that the exon/intron pattern was highly conserved between Arabidopsis and cotton, and the homeobox domain loci were also conserved between them. In addition, the expression pattern in different tissues indicated that the duplicated genes in cotton might have acquired new functions as a result of sub-functionalization or neo-functionalization. The expression pattern of WOX genes under different stress treatments showed that the different genes were induced by different stresses. CONCLUSION: In present work, WOX genes, classified into three clades, were identified in the upland cotton genome. Whole genome and segmental duplication were determined to be the two major impetuses for the expansion of gene numbers during the evolution. Moreover, the expression patterns suggested that the duplicated genes might have experienced a functional divergence. Together, these results shed light on the evolution of the WOX gene family, and would be helpful in future research.

GhCaM7-like, a calcium sensor gene, influences cotton fiber elongation and biomass production.

Calcium signaling regulates many developmental processes in plants. Calmodulin (CaM) is one of the most conserved calcium sensors and has a flexible conformation in eukaryotes. The molecular functions of CaM are unknown in cotton, which is a major source of natural fiber. In this study, a Gossypium hirsutum L.CaM7-like gene was isolated from upland cotton. Bioinformatics analysis indicated that the GhCaM7-like gene was highly conserved as compared with Arabidopsis AtCaM7. The GhCaM7-like gene showed a high expression level in elongating fibers. Expression of ß-glucuronidase was observed in trichomes on the stem, leaf and root in transgenic Arabidopsis plants of a PROGhCaM7-like:GUS fusion. Silencing of the GhCaM7-like gene resulted in decreased fiber length, but also caused reduction in stem height, leaf dimensions, seed length and 100-seed weight, in comparison with those of the control. Reduced expression of the GhCaM7-like gene caused decreased Ca2+ influx in cells of the leaf hypodermis and stem apex, and down-regulation of GhIQD1 (IQ67-domain containing protein), GhAnn2 (Annexins) and GhEXP2 (Expansin). These results indicate that the GhCaM7-like gene plays a vital role in calcium signaling pathways, and may regulate cotton fiber elongation and biomass production by affecting Ca2+ signatures and downstream signaling pathways of CaM.